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Sources of uncertainties in perturbative calculations, tadpole improvement and its role in lattice perturbation 
theory, and six recent calculations are discussed. 



1. INTRODUCTION 

Perturbation theory (PT) appears in many im- 
portant theoretical and practical roles in lattice 
field theory. Conceptual issues related to the con- 
tinuum limit can be studied in a perturbative 
framework, such as the renormalizability of lat- 
tice gauge theories 0. Perturbative calculations 
are used in relating lattice quantities to those de- 
fined in continuum schemes and in constructing 
Symanzik-improved lattice actions and operators. 
A recent determination of the strong coupling 
constant Q and a new construction of lattice chi- 
ral fermions || relied on perturbative expansions. 
PT often helps us to better understand our nu- 
merical computations and provides useful checks. 

Lattice perturbation theory is a very broad 
topic. Here, I focus on weak-coupling PT in zero- 
temperature lattice QCD. First, two important 
sources of uncertainties in lattice PT calculations 
are briefly discussed. Next, tadpole improvement 
and its role in lattice PT are described. I then 
survey six selected lattice PT papers or series of 
papers from the recent past. Developments con- 
cerning renormalons are omitted, as these will be 
presented in the talk by Sachrajda 

2. SOURCES OF UNCERTAINTIES 

In this section, I briefly discuss two important 
sources of uncertainties in any calculation which 
relies on lattice perturbation theory; namely, the 
very applicability of PT and the choice of expan- 
sion parameter in finite-order approximations. Of 
course, these issues must also be confronted in 
continuum perturbation theory. 



2.1. Applicability 

Because of asymptotic freedom, one generally 
expects that quantities in QCD dominated by mo- 
mentum scales fj, ^> Aqcd may be reliably deter- 
mined using perturbation theory. Unfortunately, 
it is very difficult to anticipate the size of nonper- 
turbative contributions in any calculation. Warn- 
ings concerning the enhancement of nonpertur- 
bative effects as a — > in the presence of power 
divergences have been issued S. The Lepage- 
Mackenzie q* scale (see below) can also serve as a 
cautionary measure of PT's reliability. Neverthe- 
less, no harm is done by trying out perturbation 
theory to see how well it works using as many 
cross-checks as possible; in fact, much might be 
learned. 

Good agreement between perturbative and 
Monte Carlo estimates of the critical quark mass 
for Wilson quarks and various Creutz ratios of 
Wilson loops was demonstrated in Ref. ||. Es- 
timates of light-quark current renormalizations 
from boosted PT M and nonperturbative chiral 
Ward identities have been recently compared 
for the improved clover fermionic action (without 
tadpole improvement). The results are summa- 
rized below (Zy, Za, Zps, and Z$ are the vector, 
axial, pseudoscalar, and scalar current renormal- 
izations, respectively). 





P 


Chiral WI 


1-loop BPT 


Zv 


6.0 


0.824(2) 


0.83 




6.2 


0.817(9) 


0.83 


Z A 


6.0 


1.09(3) 


0.97 




6.2 


1.045(10) 


0.97 


Zps/Zs 


6.0 


0.60(2) 


0.71 




6.2 


0.649(9) 


0.72 



Other nonperturbative determinations of these 



currents have also been compared to the above 
PT results in Ref. jl0| . Perturbative estimates of 
the power-divergent energy shift aEo and heavy- 
quark mass renormalization Z m in NRQCD have 
been used |yj to determine the T mass aMy^ 

from simulation results aE~r using aM^ = 
2(Z m aM — aEo) + aEy. The results compare 
very well with the nonperturbative determina- 
tions aMj^ obtained by fitting to the dispersion 

aE~r(p) = aEr +a 2 p 2 /(2aM^ '), as shown in the 

table below. 
aM aE r Z m aE aM^ a) aM^ 
1.71 0.453(1) 1.20 0.32 3.92 3.94(3) 
1.80 0.451(1) 1.18 0.31 4.08 4.09(3) 
2.00 0.444(1) 1.16 0.30 4.48 4.48(4) 

2.2. Choice of expansion parameter 

A major source of ambiguity in perturbative 
calculations is the choice of expansion parame- 
ter. Consider the perturbative prediction for an 
observable 0(Q) depending on a single external 
momentum Q of the generic form 

0(Q) = c [l + Cl (Q/ M ,RS) a s (^RS) 

+c 2 (Q/m,RS) a2(/i,RS) + ■••]. (1) 

The coefficients Cj(Q//x,RS) and the expansion 
parameter a s (/i,RS) depend on the choice of scale 
\x and the renormalization scheme (RS), i.e., the 
definition of a s (/i,RS), in such a way that the ob- 
servable 0(Q) is independent of both \x and the 
RS. However, in practice, one is forced, either by 
exhaustion or the asymptotic nature of the series, 
to truncate the perturbation expansion at some 
finite order in the coupling. Unfortunately, the 
resulting finite-order approximants are no longer 
independent of [i and the RS. 

For a good choice of expansion parameter, the 
uncalculated higher-order terms should be small. 
Since these contributions are unknown, one must 
ultimately rely on a guess about their sizes. This 
guess can be based on the sensitivity to changes 
in the scale and the RS, the apparent rate of con- 
vergence of the perturbation series for a variety 
of quantities, the size of the coefficient of the last 
term in the truncated series, or some physical cri- 
teria. The signal for a poor expansion parameter 
is the appearance of large higher-order coefficients 



in the perturbative expansions of numerous ob- 
servables. 

Choosing an expansion parameter involves fix- 
ing the RS (defining the coupling) , specifying the 
scale at which to evaluate the coupling (in the 
case of a running coupling), and determining the 
numerical value of the coupling at that scale. Var- 
ious examples of expansion parameter choices are 
described below. 

2.2.1. Coupling definitions 

Some couplings which have appeared in the lit- 
erature include: 

(a) The bare lattice coupling chq. This is a poor 
expansion parameter, yielding first-order correc- 
tions for various short-distance quantities which 
are consistently much too small in comparison to 
simulation measurements. Large expansion co- 
efficients routinely appear in perturbative series 
when expressed in terms of this coupling. 

(b) The "boosted" coupling a& = 3ao/( r Tr{7 p i a q) 
defined by rescaling the bare coupling using the 
mean plaquette. Here, the mean field contribu- 
tions responsible for the problems in (a) are ab- 
sorbed into the coupling, resulting in better be- 
haved perturbative series. 

(c) The coupling asF(q) defined through the 
Schrodinger functional using a recursive finite- 
size scaling technique |l2|,|l3| . This coupling runs 
with the finite system size L = 1/q. An alterna- 
tive coupling ctTp(q) based on the correlations of 
Polyakov loops in systems with twisted boundary 
conditions has also been proposed Jl4j . 

(d) The coupling a q q(q) = |r 2 F(r) for large 
q = 1/r defined in terms of the interquark force 
F(r) = dV/dr ||. 

(e) The coupling ay (q) defined in terms of the 
short-distance static quark potential V{q) using 
V(q) = —4:TrCFav{q)/q 2 , for q large [[§. By ab- 
sorbing higher-order contributions to the static- 
quark potential into the coupling, it is hoped 
that higher-order contributions to other physi- 
cal quantities in terms of this coupling will then 
be small. This definition facilitates scale set- 
ting since the running coupling's argument can be 
easily related to a gluon momentum. Order-ay 
agreement of perturbative with simulation results 
for several short-distance quantities has also been 



demonstrated. 

(f) The coupling a^jg- defined in the familiar 
modified minimal subtraction scheme. 



2.2.2. Scale settings 

Various scale setting prescriptions have been 
devised: 

(a) Since ceo and a& do not run, no scale setting 
is required for these couplings. 

(b) A very simple possibility is to somehow 
guess the scale. For some quantities, it should 
be possible to crudely estimate the scale likely to 
dominate the processes involved. 

(c) When the two-loop contribution is known, 
the scale may be chosen so that the one-loop co- 
efficient vanishes [fl3[. This procedure works well 
when relating couplings in different schemes (and 
is equivalent to choosing the relative scale as the 
ratio of A-parameters) but has not been tested on 
other quantities. 

(d) Another scale setting scheme is the Lepage- 
Mackenzie q* prescription ||. For a one- loop con- 
tribution Jd 4 q £(q), where q is the momentum 
of the exchanged gluon, one chooses q* such that 
a v (q*) J d 4 q £(q) = J d 4 q a v {q) £(<?)■ Inserting 
the lowest-order form for the running coupling 
into the integral on the right-hand side yields the 
result ln(g* 2 ) = Jd i q\n(q 2 ) £(q)/ Jd i q £(q). This 
procedure is the lattice analogue of the Brodsky- 
Lepage-Mackenzie prescription in continuum per- 
turbation theory[fl6||. Difficulties with this proce- 
dure can arise when J d 4 q « 0, and note that 
the mean value theorem guarantees that q* will 
satisfy < aq* < 2ir only if £(q) > for all q 
throughout the region of integration (or < for 
all such q). 



2.2.3. Value determinations 

Numerical values for the chosen coupling may 
be assigned in various ways: 

(a) In a given simulation, determining the val- 
ues for czq and otb is straightforward. 

(b) The value of the chosen coupling may be 
obtained using a perturbative expansion in terms 
of the bare couplings «o or a^. For example, in 



SU(3): 

ctv(s/a)=ao + (6. 71 — 1.751ns) «□ + • • • , 
= a 6 + (2.52-1.75his) ag + - ■ •, 
ajTg(s/a) = ao + (5. 88 — 1.751ns) a 2 , 

+ (43.41 -21.89 lns + 3. 06 ln 2 s) ag+- • • , 

= a 6 + (1.69-1.75ms) a\ 
+ (6.31-7.23 lns + 3.061n 2 s) a|+---, 

asjr(s/a)=ao + (4.62-1.751ns) a^H , 

= a h + (0.43-1.751ns) a 2 + ---. 

(c) The coupling strength can also be deter- 
mined by accurately measuring in a simulation 
some short-distance quantity whose perturbative 
expansion is reliable and known to sufficiently 
high order. In Ref. ||, the logarithm of the mean 
plaquette was used: 

-ln(±Tr[/ plaq ) = f a v (3Al/a) {l-1.185av} ,(2) 

where the scale 3.41/a is determined using the q* 
prescription mentioned above. Values of the cou- 
pling at other large values of q are then obtained 
using the familiar two-loop perturbative evolution 
equation. 

(d) A finite-size scaling procedure can be used 
to measure cusf (l) m terms of Sommer's scale r 
fTlj ] without recourse to perturbation theory. 

3. TADPOLE IMPROVEMENT 

Tadpole improvement (TI) refers to the sim- 
ple procedure of modifying any lattice gauge field 
operator by rescaling the link variable U^x) by 
a mean field factor: U^(x) — > U^(x)/uq, where 
a convenient gauge-invariant choice for the mean 
field parameter is u = (^Trf/p^) 1 / 4 §. The 
purpose of this procedure is to assist in the con- 
struction of improved lattice operators, that is, 
lattice operators with diminished discretization 
errors and lattice-to-continuum renormalization 
factors nearer to unity. 

When constructing any lattice operator 0[U] 
which is a sum of various basic gauge-field op- 
erators Oj[U], that is, 0[U] = T,j c j(a) °A U \i 
where each Oj is a simple gauge-invariant product 
of link variables, the coefficients Cj must some- 
how be determined. Lattice perturbation the- 
ory is often called upon to fix these coefficients. 



However, the Cj's generally contain large mean- 
field contributions; to reliably account for these 
tadpole effects using low-order perturbative ex- 
pansions, even when using a good expansion pa- 
rameter, is asking much from perturbation the- 
ory. Tadpole improvement offers a better alterna- 
tive: use mean-field theory instead of PT to treat 
the tadpole contributions. By tadpole improving 
the basic operators, 0[U] = X^Cj(ff) Oj[U/u Q ], 
the mean-field effects are removed nonperturba- 
tively, resulting in smaller leftover coefficients Cj 
which should be much more reliably estimated by 
low-order perturbative expansions. Note that uo 
is a nonperturbatively measured parameter when 
simulating, but its perturbative expansion must 
be used when computing the 5j in PT. 

Hence, tadpole improvement should be viewed 
as a simple means of combining mean-field theory 
with perturbation theory in order to determine 
the parameters in a lattice operator. For example, 
consider the Symanzik-improved gluon action of 
Liischer and Weisz with the chair coupling 
set to zero. The ratios of the rectangle f3 r t and 
parallelogram /3 pg couplings to the plaquette /3 P \ 
coupling change under TI as follows: 

-20^ = (1 + 2.02a s ) 2 (1 + 0.48a s ) 

P P i 

-^EI = o.03a s H> 0.03u^ 2 a s , 
P P i 

where uq is measured in the simulation (the basic 
operators without TI correction are multiplied by 
these coefficients) . This is not the same as simply 
using the boosted cib- Note that mean-field cor- 
rections are sometimes large; heavy quarkonium 
spin splittings are dramatically underestimated 
by a factor of 1/2 using tree-level couplings if TI 
is not implemented. 

4. HIGHLIGHTS FROM RECENT PAST 

Six selected lattice perturbation theory papers 
or series of papers from the recent past are sur- 
veyed in this section. For recent developments 
concerning renormalons, see Ref. 



4.1. On the viability of lattice perturbation 
theory 

Apparent discrepancies between perturbative 
and Monte Carlo estimates of various short- 
distance quantities were shown in Ref. [0] to re- 
sult from the use of the bare lattice coupling a 
as the expansion parameter. An expansion pa- 
rameter oty(q*) defined in terms of the physical 
static-quark potential was advocated. Studying 
the expectation value of the trace of a link in Lan- 
dau gauge, the critical mass for Wilson quarks, 
and various Creutz ratios of Wilson loops, the 
authors demonstrated that such discrepancies do 
not occur when a renormalizcd coupling such as 
av{q*) is used. 

Also in this paper, the nonlinear relation be- 
tween the link operator and the gauge field was 
identified as a source of large mean-field renor- 
malizations which hamper attempts to construct 
improved lattice operators. A tadpole improve- 
ment scheme, as previously discussed, was sug- 
gested to remedy this problem. 

Lastly, the onset of asymptotic or perturbative 
scaling in lattice QCD using the standard Wilson 
action was investigated. The IP-IS mass split- 
tings in charmonium and bottomonium and the 
string tension were shown to scale well for (3 val- 
ues as low as 5.7 when expressed in terms of the 
scale parameter Ay associated with ay; scaling is 
not observed for these quantities when expressed 
in terms of Ao associated with the bare lattice 
coupling. This suggests that lattice spacings used 
in current simulations are small enough for reli- 
able studies of QCD. 

4.2. in terms of ao 

In an impressive series of papers [[L9)-[g2|, 
Liischer and Weisz have recently extended to two- 
loop order the perturbative expansion of (x-^ in 
terms of the bare lattice coupling a f° r SU(AT) 
gauge theories. By matching the perturbative ex- 
pansions of corresponding correlation functions in 
the lattice and continuum theories, they found 

a Ms( s / a ) = «o + di(s) al + d 2 (s) H , 

where 

lliV 7T 

dx(s) = Ins hfciAT, 



d 2 ( S ) = d2 (s) ___ lns+ _ +fc2 + fc3iv 2 ; 

fci = 2.135730074078457(2), 
k 2 = -2.8626215972(6), 
k 3 = 1.24911585(3). 

This calculation is part of an overall strategy for 
measuring the running coupling using lattice sim- 
ulations and a nonperturbatively defined coupling 
a$F related to the Schrodingcr functional. 

The calculation exploits the background field 
technique of de Wit. In order to place their cal- 
culation on firm theoretical grounds, the authors 
first show in Ref. ]2lJ that lattice gauge theory 
with a background gauge field is renormalizable 
to all orders in perturbation theory. The proof is 
based on the BRS, background gauge, and back- 
ground shift symmetries of the lattice functional 
integral. They find that no new counterterms are 
required in addition to those already needed in 
the absence of the background field. 

An important advantage in using the back- 
ground field technique is the fact that the rela- 
tion between and «o can be extracted solely 
from the background field 2-point function: the 
3-point vertex function need not be considered. 
This dramatically reduces the number of Feyn- 
man diagrams which must be computed. Also, 
diagrams with two external legs are much sim- 
pler to evaluate than those with three external 
legs. 

The calculation is described in detail in 
Ref. and involves 4 one- loop diagrams and 31 
two-loop diagrams (7 factorizable, 7 ring, 7 tad- 
pole, 3 diamond, 3 eye, and 4 bigmac diagrams). 
All vertex factors were generated using algebraic 
manipulation programs written in Maple. The 
Feynman loop-integrals were evaluated using a 
new, innovative position-space method described 
in Ref. |pp| . This method is based on efficiently 
evaluating the free massless propagator in coor- 
dinate space using a recursion relation which ex- 
presses the propagator in terms of its values close 
to the origin. Convergence of lattice sums is ac- 
celerated using known asymptotic forms of the 
propagator. The method is even useful for evalu- 
ating Feynman diagrams with non-zero external 
momenta. 



4.3. Tadpole- improved heavy- light lattice 
operators 

In Ref. |2^| , the lattice-to-continuum renormal- 
izations of the temporal components of the point 
(A^) and point-split (Ap s ) axial currents in the 
heavy-quark effective theory (HQET) were com- 
puted to one-loop order. The standard Wilson 
fcrmionic action with massless light quarks was 
used with several values for the Wilson r param- 
eter. This was the first HQET calculation to im- 
plement tadpole improvement and to use av(q*) 
as the expansion parameter. Writing the contin- 
uum current A c in terms of the point and point- 
split lattice currents as A c (/j,) « Z(/j,a)AL(a)/2 
Z ps (fj>a)A^ (a) /(2uq), the renormalizations for 
r = 1 were found to be 

Z{fia) = l + (-1.48 + 0.318 lna/i) a y (2.18/a), 
Z ps (na) = l + (-0.76 + 0.3181na^) a y (2.13/a). 

These results may be compared to those from 
boosted PT in which one writes A c (fi) w 
yJlK hc Z{\xa)A L {a) sw y/2n bc Z ps (pa)A^ (a) , where 
k&c = uqk, c is the boosted critical hopping param- 
eter, a.b — ckq/uq, and, for r = 1, 

Z{pa) = 1 + (-1.64+ 0.318 In afi)a b , 
Z ps (^a) = 1 + (0.12 + 0.318 lna^)a b . 

Comparisons for a few values of j3 are shown in 
the table below: 



j3 /ia 


Z/2 


y/2n bc Z 


Z ps /2u \J2n hc Z^ 


5.7 1 


0.34 


0.40 


0.48 0.55 


6.0 1 


0.37 


0.41 


0.49 0.53 


6.2 1 


0.38 


0.41 


0.50 0.53 


5.7 aq* 


0.36 


0.43 


0.51 0.57 


6.0 aq* 


0.39 


0.43 


0.52 0.55 


6.2 aq* 


0.40 


0.43 


0.52 0.54 



One sees that differences between the one-loop 
results from tadpole-improved renormalized PT 
and boosted PT can be as large as 10%. 

4.4. Couplings in NRQCD 

Nonrelativistic lattice QCD (NRQCD) is an ef- 
fective field theory designed for studying hadrons 
containing heavy quarks. The NRQCD action 
includes interactions which systematically cor- 
rect for relativity and finite-lattice-spacing errors. 
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Figure 1. The one- loop corrections Sci and Sc2 
to two NRQCD couplings for /3 = 6.0 against the 
bare lattice quark mass aM. 



The couplings strengths of these interactions can 
be determined in perturbation theory. 

In Ref. pi] , the heavy-quark mass and wave 
function renormalization, energy shift, and two of 
the couplings in NRQCD were calculated to lead- 
ing order in perturbation theory as functions of 
the bare quark mass aM in lattice units. Tadpole 
improvement was implemented and results were 
given in terms of the QCD coupling av{q*) with 
scales set using the Lepage-Mackenzie q* prescrip- 
tion. The radiative corrections 8c\ and Sc2 to 
the so-called kinetic couplings, c\ and C2, where 
Cj = 1 + Scj, are shown in Fig. p]; ci is the cou- 
pling strength of the (J2k ^k) 2 interaction, and 
C2 is the coupling of the J^k term, where 
is a covariant lattice derivative. The mass renor- 
malization and energy shift were used in a de- 
termination of the 6-quark pole mass p5| . The 
complexity of the NRQCD action necessitated a 
heavy reliance on symbolic manipulations using 
Maple in carrying out these calculations. 



4.5. Structure functions 

The renormalization constants and mixing co- 
efficients for the lowest-twist lattice operators ap- 
pearing in the Wilson expansion of the product of 



two hadronic currents have recently been calcu- 
lated to one-loop order by two groups in Refs. p6| - 
p8| . These quantities are needed to extract from 
simulations the hadronic matrix elements relevant 
for determining moments of the quark and gluon 
distributions inside hadrons, i.e., the deep inelas- 
tic structure functions. The lattice operators con- 
sidered were of the form 



Q%. 

G T , .. 



= 7ti Dr 2 ■ ■ ■ Dr n X J V, 

= 2- n ^ 7ti 7 5 Dt 2 ■ ■ ■ L>r„ X f tp, 



ignoring trace terms, where "0 is a quark field, 
is a covariant lattice derivative, F^ v is the 
cloverleaf gluon field strength tensor, and A-^ are 
flavour matrices. The rank- two flavour-singlet 
operators Q M „ and G^ v mix and are related to 
the first moments of the quark and gluon distribu- 
tions, respectively; the rank-three Q^ VT is related 
to the second moment of the quark distribution. 
The quantities calculated were the renormaliza- 
tion factors 

Zu{[ia) = Ski - -r-C F (jkl In^a + B k i) (3) 

relating the bare lattice operators to a set of finite 
operators Ok{n) = Zki{^a)Oi{a) renormalized by 
requiring that the matrix elements of O in ex- 
ternal massless quark and/or gluon states having 
momentum p 2 — fi 2 are identical to the tree-level 
matrix elements of the bare lattice operators, for 
\xa <C 1 to minimize discretization effects. Only 
operators which cannot mix with lower dimen- 
sional operators due to symmetry considerations 
were studied. 

In Refs. |2f|^7|], results were obtained for both 
the Wilson and the improved clover actions in 
the chiral limit; tadpole improvement was not 
implemented. For the rank two operators, in- 
ternal quark loops were taken into account, and 
flavour singlet and non-singlet renormalizations 
were studied; only the quenched theory was used 
for the rank three quark operators. The calcula- 
tions relied heavily on the use of the algebraic ma- 
nipulation languages Form and SCHOONSCHIP, 
with modifications to properly treat the Dirac 



matrices. The following operators were consid- 



ered: 1 = Q {14} , 2 = G {14} , 3 = Q^ s 4} , 

Oa — Qui — \{Q±22 + Q433), Ob = Qui + 
Q114 — ^(Q242 + Q224 + Q343 + Q334), and Oc = 
Q{i23}j where {• • •} denotes symmetrization and 
the superscript NS denotes non-singlet (opera- 
tors are flavour singlets unless otherwise indi- 
cated). Some selected results for r = 1 (quenched) 
are given in the table below (results for the im- 
proved fermionic action are indicated by the su- 
perscript /): 



k,l 


ju 


Bu 




1,1 


16/3 


-3.165 


-15.816 


A, A 


13/3 


-18.824 


-27.389 


A,B 


2 


-0.924 


-3.603 


B,A 


4 


-2.955 


-11.803 


B,B 


19/3 


-17.540 


-18.538 


C,C 


25/3 


-19.005 


-29.815 



The rank two renormalizations using the Wilson 
action were compared with previous determina- 
tions: some discrepancies were found. 

The Wilson action results above were con- 
firmed in Ref. [E8| (except the last row which was 
not computed). However, these authors advocate 
the use of the operators Oa = Oa + Ob and Ob = 
Ob — 20 'a in order to diagonalize the anomalous 
dimension matrix. They also considered the fol- 
lowing operators: 0±b = Qa4-^(Qu + Q22 + Q33), 
O i = Q\ 1 O b = Q{ 214} , and 6 = Qf 2{1]4} , where 
[• • •] denotes antisymmetrization. The renormal- 
izations using the r = 1 Wilson action are given 
below: 



k 7fefc 

16 16/3 

4 

5 25/3 

6 7/3 



Bkk 
-1.892(06) 
15.795(03) 
-19.560(10) 
-15.680(10) 



Mathematica and Maple were used to perform 
the calculations. 

4.6. Stochastic perturbation theory 

Recently, an innovative numerical technique 
29 1 for obtaining weak-coupling perturbative ex- 
pansions of local observables in lattice QCD was 
proposed. An exciting aspect of this method is 



that it allows one to obtain much longer expan- 
sions than presently possible using conventional 
diagrammatic approaches. 

The method is based on Parisi-Wu stochastic 
quantization in which the gauge field is viewed 
as a random variable which evolves according 
to the Langevin equation. One step t — > t + e 
in the discrete Langevin equation consists of a 
sweep through the lattice, updating links using 
U u (x; t + e) = exp[— F u (x; t)]U u {x; t), where the 
driving function F n depends on the link vari- 
ables in all plaquettes containing the link between 
sites x and x + a/t, the Langevin time step e, 
and a noise matrix. One then writes J7 M (a;;t) = 
exp[A u (x;t)/y/J3\, rescales the time step e = r//3, 

expands A u {x; t) = J2 k >o P~ k/2 A ^ ( x \ f ) and ^ 
as power series in 1/ y/]3, and truncates to some 
order to transform the Langevin equation into 
a system of coupled stochastic finite-difference 
equations. The coefficients of the perturbative ex- 
pansion of any local observable W are then given 
by expectation values of composite operators of 

(k) 

the AfT' which are obtained by averaging over 
the Langevin history: 

W = EnP- n/2 (On), 
1 T 

(On) = lim =X)0«(t)- 



(4) 
(5) 



The discrete Langevin equation has O(e) system- 
atic errors, making it necessary to extrapolate the 
results to e — * 0. Stochastic gauge fixing is nec- 
essary for an acceptable signal-to-noise ratio and 
finite size errors must be controlled. 

In Ref. [p0|, the mean plaquette in SU(3), 



P = 1 



1 



(Tr U p ) 



E 

71 = 1 



(6) 



was computed to eight-loop order. Results were 
obtained on an 8 4 lattice with periodic bound- 
ary conditions for step sizes r = 0.01, 0.015, 0.02, 
then extrapolated to r — > 0. The first three co- 
efficients were found to agree with known values 
obtained analytically in Ref. |31| using the stan- 
dard diagrammatic approach, as shown in the ta- 
ble below. 





T ;p n ctpvi n 

I if li I I c_\ . v J.li. 


A rialvtir 


Cl 


1.998(1) 


2 


C2 


1.218(1) 


1.212(7) 


C3 


2.940(5) 


2.9605 


C 4 


9.28(2) 
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5. CONCLUDING REMARKS 

The evaluation of Feynman diagrams in lattice 
PT is difficult because the Feynman integrands 
are usually complicated functions of the loop and 
external momenta. Standard tools, such as Feyn- 
man parameters and partial integration methods, 
are not very helpful on the lattice. Progress 
is being made, however, with the development 
of new techniques, such as those proposed in 



Rcfs. 19|-|22|, an increasing reliance on and ex- 
pertise in using analytical computer programs, 
such as Maple and Mathematica, and the in- 
troduction of exciting new stochastic methods. 
Lattice perturbation theory continues to evolve 
and to play an important role in lattice field the- 
ory. 
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